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INTRODUCTION 


Breast  cancer,  the  most  prevalent  cancer  among  women,  has  proven  to  be  a  devastating  disease. 
The  pain  of  being  diagnosed  with  the  disease  is  compounded  by  the  lack  of  a  definitive  treatment. 
Currently,  the  most  widely  used  post-operative  hormonal  therapeutic  is  the  antiestrogen  tamoxifen 
(TAM).  TAM  acts  at  least  in  part  as  a  competitive  inhibitor  of  estrogen,  a  hormone  which  signals  cell 
growth  in  breast  cells.  Although  effective  in  postmenopausal  patients  with  estrogen  receptor-positive 
(ER+)  tumors,  its  effectiveness  in  premenopausal  patients  and  in  patients  with  estrogen  receptor¬ 
negative  (ER-)  tumors  is  substantially  lower.  Additionally,  TAM  resistance  and  cancer  recurrence 
strikes  almost  all  patients  with  advanced  breast  cancer  undergoing  TAM  therapy1.  The  serious 
limitations  of  this  current  treatment  highlight  the  need  for  a  greater  understanding  of  both  estrogen  and 
TAM  action. 

A  consequence  of  anti-cancer  therapy  is  an  alteration  of  carbon  flows  through  metabolic 
pathways.  Changes  in  metabolic  fluxes  allow  cells  to  provide  cellular  components  to  accommodate 
increased  growth.  In  human  breast  cancer  cells,  increased  growth  and  consequent  metabolic  changes 
are  initiated  by  estrogen  binding  to  the  estrogen  receptor.  Conversely,  TAM  is  a  cytostatic  drug  which 
appears  to  compete  with  endogenous  estrogen  for  binding  to  the  ER,  preventing  communication  of  the 
growth  signal.  The  effect  of  this  broken  communication  on  metabolism  is  presently  not  known.  The 
proposed  research  uses  novel  methodology  to  identify  the  specific  metabolic  pathways  and  the 
individual  metabolic  reactions  affected  by  estrogen  and  TAM  in  both  ER+  and  ER-  cells.  The 
metabolic  effects  of  these  agents,  manifested  as  changes  in  metabolic  fluxes,  will  be  quantified  with  the 
aid  of  NMR  techniques.  This  will  enable  a  description  of  the  relationship  between  primary  metabolism 
and  growth  of  cells  in  response  to  estrogen  and  TAM  treatment  to  be  developed.  This  new  information 
will  help  guide  the  development  of  new  therapies  and  new  therapeutic  agents  for  arresting  the 
metabolism  of  breast  cancer  cells  and  tumors.  The  specific  purposes  stated  in  the  original  proposal 
were  as  follows: 

•  Identify  the  primary  metabolic  pathways  that  are  operative  in  ER+  and  ER-  cells.  This  will  be 
achieved  by  in  vivo  13C  NMR  tracer  studies  and  extracellular  measurements  of  ER+  and  ER-  cells 
grown  in  a  hollow  fiber  bioreactor  (HFBR). 

•  Quantify  the  primary  and  secondary  metabolic  fluxes  in  breast  cancer  cells  in  the  presence  and 
absence  of  TAM  and  estrogen,  including  those  effects  which  may  be  reversed  by  increased 
concentrations  of  estrogen  (“ estrogen  rescue”),  and  those  effects  of  TAM  that  are  not  related  to 
ER  binding.  Flux  and  label  balance  equations  will  be  developed  for  carbon  pathways.  Using  both 
extracellular  assays  and  in  vivo  NMR  measurements  of  cells  grown  in  an  HFBR,  these  equations 
will  be  solved  to  determine  metabolic  fluxes  in  response  to  estrogen  and  TAM  treatment.  Fluxes 
through  branchpoints,  indeterminable  by  extracellular  measurements  alone,  can  be  quantified  using 
this  methodology. 

•  Examine  modes  of  estrogen  and  tamoxifen  action  by  investigating  the  relationship  between  fatty 
acid  synthesis  and  the  pentose  phosphate  pathway,  and  other  pathways  requiring  NADP+/NADPH. 
We  will  examine  the  effects  of  impaired  fatty  acid  synthesis  (resulting  from  the  addition  of  the  fatty 
acid  synthase  inhibitor,  cerulenin)  on  other  metabolic  pathways  in  the  presence  of  estrogen  and 
tamoxifen. 
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In  this  report  period,  we  have  made  significant  progress  toward  all  three  goals.  Additional 
experiments  to  determine  the  metabolic  requirements  of  ER+  MCF7  cells  have  been  performed  in 
microcarrier  culture.  Microcarrier  culture,  flask,  and  roller  bottle  studies  were  performed  in  order  to 
develop  and  demonstrate  the  ability  to  quantify  breast  cancer  cell  metabolism  under  varied  treatment 
conditions.  Metabolic  fluxes  were  calculated  for  estrogen-rescued  cells  using  NMR  data  from  cell 
extracts  and  our  metabolic  model.  We  have  determined  that  flux  analysis  using  cell  extract  NMR  data 
is  significantly  more  accurate  than  that  using  in  vivo  data;  thus  all  flux  determinations  were  based  on 
cell  extract  data. 

BODY 


A.  Methods  and  Procedures 


1.  General  Methods 

Cell  Cultures.  Human  breast  cancer  cell  lines  T47D  and  MCF7  were  obtained  from  American 
Type  Culture  Collection  and  maintained  in  Dulbecco's  Modified  Eagle's  Medium  (DMEM)  with  5% 
fetal  bovine  serum  (FBS)  at  37°C  in  a  humidified  atmosphere  with  5%  C02.  Confluent  Falcon  tissue 
culture  flasks  were  subcultured  weekly  following  trypsinization  with  0.05%  trypsin  in  saline 
/ethylaminediaminetetraacetic  acid  (EDTA).  Cell  counts  were  obtained  using  a  Coulter  Multisizer  or 
Neubauer  hemacytometer. 

Glucose,  Lactate,  and  Amino  Acid  Measurements.  Glucose  and  lactate  concentrations  were 
determined  using  a  Yellow  Springs  Instrument  Co.  2700  Select  Biochemistry  Analyzer.  Amino  acid 
concentrations  were  determined  using  an  Applied  Biosystems  130 A  high  pressure  liquid 
chromatography  system  (HPLC)  and  420 A  derivatizer  equipped  with  a  920  data  module. 

Roller  Bottle  Culture.  Cells  were  grown  in  850  cm2  Falcon  or  Coming  roller  bottles.  After 
seeding  at  15,000  -  20,000  cells/cm2,  the  bottles  were  rotated  at  1/3  RPM  for  twenty-four  hours  and  at  1 
RPM  thereafter2.  The  headspace  of  the  bottles  was  purged  with  5%  C02  in  air. 

2.  TAM  Effects  in  MCF7  Microcarrier  Cultures 

Cells  were  seeded  in  spinner  flasks  containing  0.67  g/L  Cytodex  3  microcarriers  (Pharmacia) 
and  50,000-60,000  cells/ml  (15,625-18,750  cells/cm2).  Cultures  were  initiated  in  2/5  of  the  final  culture 
volume  and  stirred  for  2  minutes  every  half  hour  for  the  first  five  hours.  After  this  time,  cultures  were 
brought  up  to  a  final  volume  of  500  ml  and  stirred  constantly  at  60  RPM.  Culture  media  was  replaced 
with  TAM  or  vehicle  control-containing  media  24  hours  after  initiation  (day  0).  Samples  were  removed 
daily  for  determination  of  cell  number.  Three  samples  were  removed  for  each  time  point.  Released 
nuclei  stained  with  crystal  violet  were  counted  on  a  hemacytometer4.  Media  samples  were  also 
removed  regularly  for  glucose,  lactate  and  amino  acid  measurement.  Cell  pellets  were  prepared  for  cell 
cycle  analysis  using  flow  cytometry.  Flow  cytometry  measurements  of  propidium  iodide-stained 
nuclei  were  obtained  using  a  Coulter  EPICS-XL  flow  cytometer. 

3.  Cell  Extraction 

Cells  were  removed  from  either  bottles  or  flasks  by  trypsinization  and  then  pelletted  by 
centrifugation  at  1200  RPM  for  5  minutes.  Cell  pellets  were  treated  with  1  ml  5%  perchloric  acid 
(PCA)  per  250  pi  of  cell  volume.  Cell  solids  were  removed  by  centrifugation  and  the  supernatent 
neutralized  with  250  pi  of  2M  KjCOj.  Precipitated  potassium  perchlorate  was  removed  and  the  extract 
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dried  by  lyophilization5.  Extracts  were  stored  at  -4°C  and  rehydrated  in  250  |il  of  D20  just  prior  to 
NMR  analysis. 

Using  extraction  to  characterize  the  intracellular  isotope  distribution  and  subsequently 
determine  intracellular  fluxes  has  been  questioned  because  of  concerns  that  extraction  affects  the 
chemical  make-up  of  the  isotope-containing  metabolites.  However,  preliminary  measurements 
indicated  that  the  major  metabolites  observable  in  breast  cancer  cells  by  13C-NMR  are  the  amino  acids 
glutamate  and  aspartate,  both  of  which  are  unaffected  by  acid  extraction. 

There  are  many  advantages  to  the  extraction  NMR  method  when  compared  to  the  proposed  in 
vivo  methods.  High  concentrations  of  intracellular  aspartate  and  glutamate  can  be  attained  that  produce 
very  informative  spectra  on  a  high  powered  500  MHz  spectrometer.  In  these  spectra,  isotopomers,  or 
multiply  labeled  species,  can  be  discerned.  A  measured  isotopomer  distribution  can  be  used  to  increase 
the  accuracy  of  the  calculated  fluxes.  In  addition,  numerous  flasks  or  bottles  can  be  grown 
simultaneously  prior  to  extraction.  Investigating  different  treatments  and  controls  concurrently  allows 
efficient  comparison  of  treatments  and  verification  of  their  metabolic  consequences. 

4.  Metabolic  Response  to  Estrogen  Rescue  Treatment  Using  1-13C  NMR  and  Cell  Extraction 

Six  roller  bottles  were  seeded  with  MCF7  cells.  Samples  of  the  bulk  media  were  taken  every  8 
hours  to  determine  the  concentrations  of  glucose,  lactate,  and  the  amino  acids.  The  experiment  lasted 
for  8  days.  The  bottles  were  fed  medium  containing  5  pM  tamoxifen  on  the  first  and  third  days.  On 
the  fifth  day  media  containing  1-13C  glucose  was  supplied  to  all  six  bottles.  Three  bottles  also  received 
0.5  pM  estrogen.  This  cell  cycle  hold-up  followed  by  stimulation  is  the  basis  of  estrogen  rescue.  On 
the  seventh  day  the  cytoplasm  of  all  the  cells  in  each  of  the  six  bottles  was  extracted.  All  six  bottles 
were  trypsinized  and  the  cell  pellet  resuspended  in  8  ml  of  conditioned  medium.  Conditioned  medium 
was  used  to  limit  transport  of  unlabeled  metabolites  into  the  cells.  A  small  portion  of  this  resuspension 
(0.2  ml)  was  used  for  cell  density  determination  on  a  Coulter  Counter.  The  remainder  of  the  cells  were 
re-centrifuged  and  used  for  extraction. 

5.  NMR  Analysis 

13C  spectra  were  acquired  on  a  500MHz  Bruker  DRX  spectrometer  with  a  5  mm  broad-band 
probe.  A  ninety-degree  pulse  was  used  with  a  6  second  interscan  delay.  Spectra  were  acquired 
overnight,  resulting  in  the  signal  averaging  of  approximately  5600  scans.  Concentrations  of  labeled 
species  were  determined  by  calibration  to  a  2  mM  ethylene  glycol  peak. 

6.  Development  of  a  Biochemical  Pathway  Model  and  Mathematical  Description 

Measured  external  metabolite  fluxes  and  fractional  enrichments  were  inputted  into  a 
mathematical  description  of  biochemical  pathways  to  determine  the  unknown  internal  fluxes.  In  batch 
growth  (flask  and  roller  bottle  culture),  external  fluxes,  fiex  (pmol/cell-hr),  are  calculated  from  linear 
time-course  concentration  profiles  using  equation  l6: 

f  =l.^L=h. 

1,ex  N  dt  N  (!) 

where  C,  is  the  concentration  of  metabolite  i,  N  is  number  of  cells  and  h  is  the  slope  of  the  linear 
concentration  profile.  Currently  the  metabolically  important  substrates  that  are  measured  are  glucose, 
lactate,  glutamate,  glutamine,  aspartate,  serine,  and  alanine.  Fractional  enrichments  are  determined  by 
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integrating  decoupled  carbon  spectra  to  determine  the  concentration  of  label,  Liy,  on  atom  y  in 
metabolite  i.  The  total  metabolite  concentration  in  the  extract,  C,  was  measured  by  HPLC  or  enzyme 
assay  and  was  combined  with  the  concentration  of  label  to  calculate  the  absolute  fractional  enrichment, 


(2) 


Conversion  of  fractional  enrichments  and  external  flux  data  into  internal  fluxes  requires  the 
development  of  a  biochemical  pathway  model.  It  is  impossible  to  include  all  of  the  many  pathways  in 
a  single  model.  Pathways  retained  in  the  model  have  either  significant  flows  of  carbon  or  considerable 
effects  on  the  isotope  distribution.  The  current  model  is  shown  in  Figure  1  (see  next  page)  and 
contains  the  following  pathways:  glycolysis,  the  TCA  cycle,  the  pentose  phosphate  pathway, 
glutaminolysis,  malate  transport,  tricarboxylate  transport,  malic  enzyme,  and  pyruvate  carboxylase. 

Fluxes  are  calculated  using  metabolite  and  isotope  mass  balances.  Because  metabolic  steady 
state  is  assumed7,  the  intracellular  accumulation  is  negligible,  and  metabolite  balances  simply  equate 
the  sum  of  the  fluxes  producing  a  metabolite  to  the  sum  of  the  fluxes  consuming  it.  All  of  the  balances 
derived  from  Figure  1  are  given  in  Table  1: 


Table  1  Metabolite  balances  derived  from  Figure  1 


f|  =  f2  +  f3 

fn  +  4  =  f,4 

f,  +  2/3  f,  =  f4 

f.4  =  f,3  +  f16 

f4  =  2  f5 

fn  +  fn  +  f22  =  fi7  +  ^8  +  f2i  +  f25 

f,  +  1/3  f,  =  f6 

f|7  +  f)8  +  f|Q  =  f9  +  f]4 

_ £jJk±  =  fin  +  fn  +  f.9  +  f,Q _ 

f0+f„  =  f,n 

Isotope  balances  equate  the  sum  of  the  fluxes  of  isotope  into  a  single  atom  of  a  metabolite  to  the  sum 
of  the  fluxes  out. 


Figure  2.  Example  of  an  isotope  balance.  Each  circle  represents  carbon  atoms  in  one  of  three  example  molecules: 
A,  B  and  C.  The  darkened  circles  are  partially  enriched  by  13C.  The  three  carbon  atoms  of  each  molecule  are 
numbered  1-3  starting  at  the  top.  This  particular  set  of  reactions  represents  the  recombination  of  two  pathways. 

The  isotope  balance  around  the  first  carbon  of  metabolite  C  in  Figure  2  is  written  as 


fj-Al  +  f2Bl  =  f3Cl  (3) 

Where  Al,  Bl,  and  Cl  are  the  fractional  enrichments  of  the  first  carbons  of  metabolite  A,  B,  and  C, 
respectively.  The  product  of  f-C  thus  has  units  of  |imol/cell-hr.  Two  more  balances  can  be  written  for 
the  example,  one  for  the  second  carbon  and  one  for  the  third.  An  isotope  balance  in  this  form  can  be 
written  around  every  atom  of  every  metabolite  in  the  pathway  model.  In  Figure  1  there  are  44  such 
balances.  Because  very  few  of  the  enrichment  fractions  were  actually  measured,  it  was  necessary  to 


8 


G 1 -G2-G3-G4-G5-G6 


Cytoplasm 


Ru  1  -Ru2-Ru3-Ru4-Ru5-P 


9 


Clark,  Douglas  S 


reduce  the  number  of  equations.  It  was  possible  to  produce  5  equations  for  each  of  the  five  isotope 
measurements  made. 

The  method  of  least  squares  was  used  to  find  the  set  of  unknown  fluxes  by  minimization  of  the 
sum-of-squares,  F,  which  is  given  by 


r  f[y,-yiWf 


(14) 


where  y(  is  any  measured  value,  either  fractional  enrichment  or  external  flux,  and  y^x)  is  the  adjusted 
value  of  that  variable  constrained  by  the  metabolite  and  isotope  balances.  There  are  n  measurements 
and  x  represents  the  adjusted  unknown  fluxes.  The  standard  deviation  of  each  measurement,  <jr  is 
included  to  weigh  the  sum  in  favor  of  more  accurate  measurements. 

Errors  in  the  values  of  the  internal  fluxes  were  calculated  by  simulation  around  the  errors  of  the 
input  measurements8.  This  is  done  by  first  generating  a  statistically  probable  set  of  input  data  inside  a 
Gaussian  distribution  of  the  measurements.  A  new  set  of  internal  fluxes  was  calculated  from  this  new 
set  of  input  “data.”  This  process  is  repeated  M  times,  after  which  the  error  of  each  internal  flux  is 
given  by 

mJ!mS 

V  j=i  /  (15) 

A 

where  f  is  the  set  of  best-fit  fluxes,  determined  from  the  initial  values  of  the  measurements.  These 
errors  of  the  internal  fluxes  have  three  purposes:  they  establish  the  accuracy  of  the  results,  they 
discriminate  plausible  from  implausible  pathway  models,  and  they  establish  the  statistical  significance 
between  different  test  conditions. 


7.  Isotopomers  and  Simulated  Annealing 

Advances  in  experimental  techniques  enabled  the  detection  of  isotopomers  in  internal 
metabolites  by  NMR  as  shown  in  Figure  3.  That  the  concentration  of  these  multiply  labeled 
metabolites  can  be  used  in  the  analysis  to  increase  the  accuracy  of  flux  results  and  add  more  pathways 
to  the  biochemistry  model  is  a  significant  finding  of  our  research.  The  only  pathway  that  joins  together 
metabolites  in  Figure  1,  and  hence  is  able  to  form  isotopomers,  is  the  citrate  synthase  pathway. 
Because  citrate  synthase  is  contained  within  the  mitochondria  and  isotopomers  accumulate  in  extra- 
mitochondrial  pools,  additional  pathways  in  and  out  of  the  mitochondria  were  added. 
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Figure  3.  Isotopomers  of  aspartate  visible  in  a  ,3C-NMR  spectrum  of  an 
MCF7  cell  extraction  bv  Derchloric  acid  following  estrogen  stimulation 


In  order  to  use  isotopomer  measurements  to  calculate  fluxes,  the  analysis  method  had  to  be 
enhanced.  The  problem  encountered  was  that  the  desired  results,  determining  all  of  the  fluxes  in 
Figure  1,  could  not  be  calculated  from  the  available  data,  a  set  of  extra-cellular  fluxes  and  a  set  of 
internal  metabolite  isotopomer  concentrations.  There  were  two  major  limitations  in  achieving  this 
goal.  The  first  is  the  scope  of  the  problem,  i.e.  the  large  number  of  variables  and  equations.  Any 
method  for  determining  a  set  of  fluxes  is  based  on  isotopomer  balances,  which  are  similar  in  form  to 
the  isotope  balances  in  equation  3.  However,  if  a  balance  for  every  isotopomer  of  each  metabolite  was 
considered  (there  are  over  1000),  the  problem  would  be  intractable.  The  second  restriction  was  that  it 
was  not  possible  to  recast  the  isotopomer  balances  into  a  simple  analytical  form  that  would  express 
each  flux  as  a  function  of  the  isotopomer  measurements  alone: 

flux;  =  f(all  isotopomers)  (4) 

This  is  due  to  the  non-linear  relation  between  isotopomer  measurement  and  flux,  and  the  self- 
dependent  nature  of  cyclic  biological  pathways.  Therefore,  a  numerical  approach  based  on  a 
simulated  annealing  algorithm  was  used  to  address  these  issues. 

The  premise  of  the  simulated  annealing  algorithm  involved  guessing  an  initial  set  of  fluxes  and 
calculating  the  accompanying  isotopomer  pattern.  These  calculated  isotopomers  were  compared  to  the 
measured  values  and  then  a  new  set  of  fluxes  was  guessed.  This  cycle  was  repeated  until  the 
difference  between  the  guessed  and  measured  isotopomer  concentrations  reached  a  minimum.  This 
method  allows  a  global  minimum  to  be  approached  in  a  highly-dimensional  non-linear  system. 

The  intractability  of  the  system  was  encountered  when  the  isotopomer  pattern  was  calculated 
from  the  guessed  set  of  fluxes.  Without  reduction,  an  unwieldy  matrix  with  elements  equal  to  the 
square  of  all  of  isotopomers  in  the  cell  would  have  to  be  inverted.  However,  there  is  another  piece  of 
information  that  can  render  this  problem  tenable:  only  a  single  source  of  isotope  (1-13C  glucose)  was 
supplied  to  the  cells.  The  amount  of  isotope  in  any  isotopomer  is  therefore  dependent  only  on  the 
fluxes  (pathways)  that  supply  it  with  isotope.  This  is  the  basis  for  a  tracing  algorithm  that  generates  an 
isotopomer  matrix  with  only  as  many  elements  as  the  square  of  all  isotopomers  observed.  This  is  more 
that  two  orders  of  magnitude  smaller  than  the  “brute  force”  approach  and  proportionally  that  much 
faster  to  solve. 
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factor.  This  factor  scales  from  0,  irreversible,  to  1,  rapidly  reversible,  and  can  be  incorporated 
seamlessly  into  the  isotopomer  matrix. 

All  of  these  advances  in  methodology  enabled  the  analysis  of  the  isotopomer  data  and  an 
expansion  of  the  metabolic  model. 

B.  Results  and  Discussion 

1.  TAM  effects  in  MCF7  microcarrier  cultures 

MCF7  cells  were  subjected  to  1  |lM  TAM.  The  cell  number  on  day  5  for  TAM-treated  flask 
was  approximately  30%  of  the  vehicle  control  flask.  The  growth  curves  are  shown  in  Figure  4. 
Profiles  of  glucose,  lactate,  and  glutamine  concentrations  as  a  function  of  time  were  determined  and  are 
shown  in  Figure  5.  These  profiles  were  then  fit  to  a  standard  quadratic  equation.  The  derivative  of  the 
fit  was  used  to  calculate  uptake  and  production  rates  of  glucose,  lactate,  and  glutamine  respectively. 
The  results  are  shown  in  Table  2.  Cell  cycle  phase  was  determined  on  day  5  and  is  shown  in  Figure  6. 
Cell  morphology  was  documented  using  phase-contrast  light  microscopy  and  is  shown  in  Figure  7. 


Glucose 

Glutamine 

Lactate 

Control 

4.2+0.8 

3.9±0.7 

1. 4+0.3 

1  pM  TAM 

0.0±0.1 

2.5+0. 1 

4.8±0.1 

Table  2.  Uptake  rates  (mmol/day/109  cells)  for  MCF7  microcarrier  cultures,  day  5. 

The  above  experiment  indicated  that  certain  metabolic  effects  of  TAM  can  be  quantified  using 
microcarrier  cultures.  Glucose  uptake  on  a  per  cell  basis  were  significantly  lower  for  TAM-treated 
cultures  than  for  control  cultures  of  MCF7  cells.  Interestingly,  lactate  was  consumed  after  several  days 
in  culture.  This  consumption  was  higher  for  TAM-treated  cells.  Lactate  consumption  did  not  occur  in 
T47D  cultures,  for  which  data  were  presented  in  the  last  report.  T47D  cells  showed  less  sensitivity  to 
TAM  and  higher  glycolytic  capacity  than  MCF7  cells.  This  microcarrier  data  will  be  used  for 
comparison  to  flux  results  obtained  from  cell  extract  experiments.  Cell  cycle  analysis  was  used  to 
investigate  the  nature  of  TAM  action.  The  results  indicate  that  TAM  was  acting  as  a  cytostatic  agent, 
causing  cell  cycle  arrest  in  G0/G1. 


2.  Intracellular  fluxes  of  estrogen-rescued  MCF7  cells 

The  time  response  of  glucose  and  lactate  following  estrogen-rescue,  as  described  in  the 
methods,  is  shown  in  Figure  8.  The  specific  uptake  of  glucose  and  production  of  lactate  both  increased 
in  the  presence  of  estrogen.  The  NMR  spectrum  of  the  extract  is  shown  in  Figure  9.  The  areas  of  the 
labeled  peaks  were  integrated  and  inputted  into  the  algorithm  to  calculate  the  fluxes  shown  in  figure 
10. 

There  were  a  few  fluxes  that  were  significantly  affected  by  estrogen.  There  was  a  six-fold 
increase  in  the  pentose  phosphate  pathway  (PPP)  flux.  This  supports  the  hypothesis  that  estradiol 
stimulates  the  PPP  in  order  to  produce  more  NADPH,  which  is  essential  for  biosynthesis.  Pyruvate 
carboxylase  and  glutaminolysis  both  increased.  These  anaplerotic  pathways  supply  the  mitochondria 
with  TCA  cycle  intermediates.  These  intermediates  are  the  building  blocks  for  many  biosynthetic 
pathways.  Again  this  supports  the  hypothesis  that  estrogen  stimulates  the  biosynthetic  pathways.  The 
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preference  for  one  of  the  two  orientations  of  carbon  scrambling  in  the  TCA  cycle  supports  the 
hypothesis  that  all  of  the  TCA  cycle  enzymes  are  grouped  into  a  complex  that  passes  metabolites 
without  releasing  them  freely  into  the  bulk. 

3.  Preliminary  Inhibitor  Results 

Six  flasks  were  grown  in  different  concentrations  of  cerulenin  to  determine  which  concentration 
would  most  effectively  attenuate  fatty  acid  synthase  activity  without  causing  cell  death.  The  results  are 
shown  in  Figure  1 1.  A  final  working  concentration  of  26.9  |iM  was  chosen. 

C.  Recommendations 


The  Statement  of  Work  in  the  original  proposal  is  included  below.  Recommendations  in 
relation  to  the  Statement  of  Work  are  given  in  Italics. 

Task  1:  Months  1-2:  Serially  wean  MCF7  ER+,  MCF7  ER-,  and  T47D  ER+  cells  to  2%  serum. 
Concerns  were  raised  about  the  use  of  low  levels  of  serum.  MCF7  and  T47D  cell  lines  are  being 
maintained  in  5%  serum. 

Task  2:  Months  2-5:  Implement  more  precise  analytical  control  and  monitoring  capabilities  of  existing 
hollow  fiber  bioreactor  (HFBR)  support  apparatus.  A  more  robust  temperature  control  algorithm 
implementing  cascade  control  has  been  created.  This  controller  is  capable  of  maintaining  the  desired 
set  point  (37  °C)  plus  or  minus  0.1° C.  However,  this  system  will  not  be  used  since  data  obtained  from 
batch  extraction  experiments  afford  significantly  higher  accuracy.  In  an  extraction  experiment  it  is 
possible  to  obtain  a  100-fold  increase  in  NMR  sensitivity.  An  experiment  takes  one  week  as  opposed  to 
three  months  and  can  test  more  than  6  conditions  as  opposed  to  1  in  an  HFBR. 

Task  3:  Months  5-6:  Treat  HFBR  with  attachment  factors  and  test  cell  growth.  T47D  cells  were 
grown  in  an  HFBR  for  3  months  at  which  point  the  cells  where  removed  by  trypsinization.  Both  NMR 
and  culturing  techniques  indicated  that  there  was  little  viability  in  the  5g  cell  pellet.  A  modified  HFBR 
would  be  needed  to  attain  a  higher  viable  cell  density.  To  overcome  the  difficulties  with  the  HFBR 
apparatus,  roller  bottles  are  used  to  obtain  the  necessary  inputs  for  flux  analysis. 

Task  4:  Months  7-8:  Perform  growth  studies  to  determine  cerulenin  concentration  to  be  used  in 
HFBR  experiments.  This  task  has  been  completed. 

Task  5:  Month  9:  Inject  treated  HFBR  with  ER+  MCF7  cells  and  allow  them  to  grow  to  high  density. 
See  Task  3. 

Task  6:  Months  10-14:  Execute  serial  step  changes  in  medium  composition,  such  as  TAM 
concentration  and  13C-labeled  metabolites,  while  continually  measuring  internal  and  external 
metabolites.  Serial  step  changes  were  not  performed,  but  equivalent  information  was  obtained. 

Task  7:  Month  15-16:  Perform  cell  extract  experiments  to  confirm  spectral  assignment  of 
metabolites.  As  was  described  in  Materials  and  Methods,  the  extraction  technique  proved  to  be  more 
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beneficial  than  expected.  The  ability  to  detect  isotopic  distribution  by  13C-NMR  measurements  of 
extracts  prompted  us  to  use  cell  extract  methodology  for  all  subsequent  experiments. 

Task  8:  Months  17-18:  Perform  data  analysis.  Check  for  inconsistencies  with  existing  metabolic 
model.  It  was  possible  to  perform  this  task  with  13C-NMR  extract  measurements.  As  described  in  the 
results  section,  many  models  have  been  tested  and  a  simulation  method  has  been  developed  to  discern 
between  models. 

Task  9:  Month  19:  Inject  treated  HFBR  with  ER+  T47D  cells  and  allow  to  grow  to  high  density. 

See  task  3. 

Task  10:  Months  20-28:  Repeat  tasks  6-8  but  for  the  T47D  cell  line. 

Task  11:  Month  28:  Inject  treated  HFBR  with  ER-  MCF7  cells  and  allow  to  grow  to  high  density. 

See  task  3. 

Task  12:  Months  29-36:  Repeat  tasks  6-8  but  for  the  ER-  cell  line. 


KEY  RESEARCH  ACCOMPLISHMENTS 

•  Identified  key  operative  metabolic  pathways  in  ER+  MCF7  breast  cancer  cells 

•  Developed  novel  methodology  employing  NMR  data,  including  isotopomers,  from  cell  extracts  and 
a  metabolic  model  to  calculate  intracellular  fluxes 

•  Used  this  methodology  to  characterize  metabolic  responses  of  MCF7  cells  to  estrogen  and 
tamoxifen 

REPORTABLE  OUTCOMES 


Forbes,  N.S.,  Blanch,  H.W.  and  Clark,  D.S.  “Metabolic  Flux  Analysis  to  Characterize  the  Response  of 
Human  Breast  Cancer  Cells  to  Hormones,”  poster  presented  at  the  Metabolic  Engineering  II 
Meeting,  Elman,  Germany,  October  25  1998. 

Forbes,  N.S.,  Blanch,  H.W.  and  Clark,  D.S.  “Analysis  of  Metabolic  Fluxes  in  Mammalian  Cells”.  In: 
Bioreaction  Engineering:  Modelling  and  Control.  Eds.  Schugerl,  K.,  Bellgardt,  K.  Springer- 
Verlag,  Berlin,  Germany,  in  press. 

Wicklund,  K.A.,  Blanch,  H.W.  and  Clark,  D.S.  “Characterization  of  Human  Breast  Cancer  Cells  in 
Microcarrier  Culture:  Effects  of  Glucose  and  Oxygen,”  poster  presented  at  the  Whitaker 
Foundation  Annual  Meeting,  La  Jolla,  CA,  August  18, 1999. 
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Wicklund,  K.  A.,  Blanch,  H.W.  and  Clark,  D.S.  “Characerization  of  Human  Breast  Cancer  Cells  in 
Microcarrier  Culture,”  presented  at  the  American  Chemical  Society  National  Meeting,  Anaheim, 
CA,  March  25, 1999. 

Wicklund,  K.A.,  Blanch,  H.W.  and  Clark,  D.S.  “Effects  of  Glucose  and  Oxygen  on  MCF7 
Metabolism  and  Response  to  Tamoxifen.”  Manuscript  in  preparation. 


CONCLUSIONS 


In  summary,  we  have  developed  advanced  methodology  that  has  allowed  us  to  investigate  the 
metabolic  consequences  of  estrogen  and  tamoxifen  treatment.  We  have  identified  for  the  first  time 
isotopomers  following  1-13C  glucose  treatment,  which  was  enabled  by  advances  in  the  proposed 
technique.  We  have  also  developed  an  algorithm  to  analyze  isotopomer  data  that  is  fast,  general,  and 
increases  the  accuracy  and  extent  of  flux  calculations.  These  improvements  in  methodology  have 
allowed  us  to  quantify  more  completely  internal  metabolic  fluxes  in  ER+  human  breast  cancer  cells 
following  estrogen  rescue.  Results  from  this  analysis  have  confirmed  our  hypothesis  about  the  action 
of  estrogen.  Specifically,  increases  in  the  pentose  phosphate  pathway,  glutaminolysis,  and  flux 
through  pyruvate  carboxylase  were  observed  following  estrogen  stimulation.  Finally,  we  have 
completed  preliminary  investigations  using  the  fatty  acid  synthase  inhibitor,  cerulenin. 
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Figure  6:  Cell  cycle  analysis  of  MCF7  microcarrier  cultures  on  day  5 


Figure  7:  Phase-contrast  light  microscopy  image  of  MCF7 
microcarrier  control  cultures  on  day  5. 
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Figure  8.  Glucose  and  lactate  response  following  estradiol  treatment  of  MCF7  cells. 
External  fluxes  calculated  by  normalizing  the  first  derivative  of  the  concentration  vs.  time 

profile  by  the  cell  number. 
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Figure  11:  Effect  of  different  concentrations  of  cerulenin  on  the 
glucose  uptake  rate  of  MCF7  cells. 
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